Commensurability Transition and Stripe Phases in the Ginzburg-Landau Theory 
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We phenomenologically describe the thermodynamics of charge and spin density waves in doped 
high-Tc oxides. We have explicitly calculated stable non-homogeneous solutions in the incompress- 
ible spin driven stripe phase, where stripes are static soliton-like charge density waves (CDW), 
and in the vicinity of their critical point, where CDWs become harmonic. Our phase diagram 
points to a commensurability transition separating the low (LI) and high (HI) incommensurable 
phases. Besides, we demonstrate by rigorous group symmetry arguments that the stripe criticality 
is compatible with a second order phase transition. 
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I. INTRODUCTION 

Stripe formation in the high-T c oxides has been widely 
considered by many theorists in the recent years. The 
present experimental evidences point to a curious elec- 
tronic order, where carriers introduced by doping (holes) 
segregate into 'rivers of charge'. Those structures, gener- 
ically called stripes, are considered as 'nearly one dimen- 
sional' objects immersed in an insulating antiferromag- 
netic sea of spins. This picture was inferred more ac- 
curately after recent neutron diffraction experiments^, 
which observed incommensurable antiferromagnetic do- 
mains separated by charge stripes localized in the Neel 
wall regions. It is found that the magnetic scattering 
peaks define a superzone which is half the size of the one 
associated with the charge order, indicating that stripes 
introduce a 7r-phase shift between neighboring magnetic 
domains. The above structure results from competing 
effects, which in conjunction, reduce the energy cost for 
breaking antiferromagnetic bonds and minimize the ki- 
netic energy of holes. 

The incommensurability of the peaks with the recipro- 
cal lattice is a strong argument in favor of the idea that 
holes condense in stripes of finite width oriented along 
a given direction in the insulating planes. The temper- 
ature dependence of the peak splitting, showing that in 
some systems the charge-ordered peaks appear at higher 
temperatures than the magnetic onesi, also discards the 
hypothesis of nested Fermi surfaces. Currently, the stripe 
picture is considered as a good candidate for explaining 
many of the properties (normal and superconducting) of 
the high-r c and other layered oxides^. In contrast, there 
is still a lot of controversy about the nature of the cor- 
rect spin order at the Neel walls, i.e spiral or collinear. 
Both of them are possible under the available experimen- 
tal datai^, despite theoretical difficulties to justify spiral 
order for the ground stated. 

Another key point concerning the stability of stripes 
is their robustness in the absence of long-range anti- 
ferromagnetism, at temperatures considerably above of 
the three-dimensional Neel ordering. They even sur- 



vive the absence of two-dimensional antiferromagnetic 
correlations, which occurs at higher temperature, show- 
ing that the transition is charge-driven^. This phe- 
nomenon is ascribed to the separation of charge and spin 
fluctuations^, which establishes strong dynamic correla- 
tions that lower the large energy due to the hole repulsion 
inside the stripes. This behavior is quite ubiquitous in 
all oxide families, except for the chromium oxides, where 
the stripe and spin orders collapse together—. 

The current idea among theorists (originally suggested 
by Anderson^ ) is that a sufficiently large repulsive on-site 
U potential might allow for antiferromagnetism in the 
oxygen planes. The antiparallel alignment of the antifer- 
romagnetic ordering at half filling, lowers perturbatively 
the kinetic energy in doped samples, because carriers 
(holes in most cases) can hop to neighboring sites. An- 
tiferromagnetism also produce an attractive force among 
holes, in order to minimize the number of broken antifer- 
romagnetic bonds. The competition between this short 
range attraction and the long range Coulomb repulsion 
of holes, may stabilize the non-homogeneous phase of 
stripes, frustrating the tendency to macroscopic phase 
separation. From the view point of continuous models, 
the hole repulsion competes with the "helicity" of the spi- 
ral spin order , meaning that the antiphase of domains 
is crucial to achieve the thermodynamic stability of the 
stripe regime^. The other possibility, magnetic domains 
in phase, should be unstable in relation to the homoge- 
neous state. This fact, also rules out the Ising limit, with 
very narrow walls between domains. 

The physics involved in the present problem, embraces 
as much the long range scale of charge interactions as 
the microscopic scale of spin fluctuations, where theoret- 
ical considerations in the basis of a macroscopic mean 
field may seem meaningless at first sight. However, 
for the stripe regime, spin-spin and charge-charge cor- 
relation lengths should be comparable in a certain re- 
gion of doping and temperature, with non- vanishing two- 
dimensional antiferromagnetic correlations. In spite that 
long range antiferromagnetic order is not attained in two 
dimensions, typical correlation lengths are of the order of 
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100 — 200 lattice parameters'^. We then define a two- 
dimensional Tjv as the temperature where the magnetic 
neutron diffraction peaks are not resolved any moreii. 
In a similar way, a critical point T$ for the stripe or- 
der is assumed^. The above facts justify a theoretical 
treatment using the Ginzburg-Landau (GL) approach to 
elucidate the internal structure of stripes. 

The central idea of this paper is to assess the proper- 
ties of the phase diagram, shown in a qualitative fashion 
in Fig. n an d to theorize about the nature of the cor- 
responding phase transitions. We will concentrate our 
study in two regions, where controversial statements are 
found in the literature, namely region LI, for the incom- 
pressible stripe phase below Tn, and region HI, for the 
paramagnetic phase immediately below the stripe criti- 
cal point Ts- We will show that the charge distribution 
suffers a transition of commensurability, which evolves 
from the low incommensurable regime in phase LI, char- 
acterized by a typical soliton-likc distribution, to a highly 
incommensurable one, which extends until phase HI and 
collapses at T§, following possibly a second order transi- 
tion. 
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FIG. 1: Cuprate phase diagram temperature vs doping. AF 
is the true long range antiferromagnetic order, below x ~ 0.02. 
Inside the quarter of dome is the 2D (quasi-long range) an- 
tiferromagnetic phase, coexisting with the superconducting 
one (SC). LI is the spin driven Low Incommensurable phase 
of stripes. Regions I and II are High Incommensurable phases 
with and without spin order, respectively. The harmonic High 
Incommensurable phase is denoted by HI, right below the 
High Symmetry (HS) phase, where stripes disappear. To in- 
dicates a second order transition separating high and low in- 
commensurable phases. 

In the first step of the calculation, we study the 
marginal stripe phase transition below T/v, building a 
GL free energy with the staggered magnetization as the 
main order parameter (OP) and the charge density play- 
ing the role of a secondary OP. This construction is based 
on a microscopic picture, which is supplemented by sym- 
metry considerations. In contrast, in the second part of 
the paper, we study the highly incommensurable phase 
in the absence of antiferromagnetism, the transition be- 
ing driven by the charge density. We develop theoretical 
arguments entirely based on symmetry to build the GL 



free energy for the HI phase, where the CDW is har- 
monic in the vicinity of Tg. We calculate the amplitude 
and wavelength of the CDW and estimate the limits of 
the harmonic approximation. This one and other conjec- 
tured phases are displayed in the diagram of Fig.l. We 
end the paper with a discussion of our results. 



II. LOW INCOMMENSURABLE STRIPE 
PHASE 

Attention has been paid by some theorists to the ten- 
dency of attraction between domain walls at large dis- 
tance. Priadko and coworkers^ believe that this fact leads 
to an ubiquitous tendency to phase segregation that ren- 
ders the stripes unstable in the small doping limit. They 
demonstrated the existence of an attractive force between 
asymptotically separated charged domain walls when the 
charge effect is considered in first order by means of a 
chemical potential. Here, we argue that the spin ex- 
change "traps" the holes inside domain walls, leading to 
a direct competition between the spin order attraction 
and the long range repulsion of holes, which stabilizes 
the inter-stripe separation. In spite that charges play 
the role of a secondary OP, their interactions must be 
included at least in second order in the free energy. In 
the cuprate case, it is experimentally known that for the 
doping range 0.05 < x < g, the effect of doping increas- 
ing is counterbalanced by the reduction of the domain 
width d, and vice-versa, following the empirical relation 
x ~ 2- I n this phase, the stripes are internally charge 
incompressible and the low doping regime may be sim- 
ply thought as a re-scaling of the non-interacting stripe 
system. 

First of all, we propose a constraint in the form of a 
global charge conservation 

Jvrp(r) = Q , 

where the integration is performed over a two dimen- 
sional system with linear dimensions (l,l z ), with the pr- 
axis along the charge modulation direction. From now 
on, the charge of holes, the electron spin, and the lat- 
tice constant are normalized to unity and will be omit- 
ted. If the stripes are non-interacting, the system can 
be rescaled along the direction z in terms of the domain 
width d, considering that two neighboring stripes are in- 
finitely distant (i.e. d — > oo). This way, the charge con- 
servation principle can be simply restated as 

fd/2 

n = lim / dz p(z) = h{T) , (1) 

J-d/2 

where h(T) is a scalar function of temperature, h(T) be- 
ing proportional to a linear charge density Q/l divided 
by the number of domains l z /d along the modulation di- 
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rection, i.e. 



h(T) = jj-d = x d 



or i a 2- Therefore, the constraint is quite appro- 
priate for the low doping regime. 

The total energy of each plane must take into account 
the long range potential for holes along with the exchange 
interaction of spins. The simplest static Hamiltonian for 
a square lattice of N sites is 



H 



H c + H s = V 

i,5 
N 

-^2 Ji n l (a 1 Sf l + a 2 Sl i ) + 

i 

N 

^ Ji,5nin l+ s(/3iSiiSii + 8 + (3 2 S2iS2i+s) (2) 



i,5 



with rii = 0,1, the spin occupation number. The in- 
teractions are extended to distant neighbors, which are 
labeled by 8. The first term in (J5J) yields the repulsive 
interactions of (spinless) holes with coupling (I /eg). The 
second term includes the crystal field contribution, with 
anisotropy parameters a\ and ot 2 . The last term consid- 
ers the exchange interactions among spins (Ji,s > 0), 
which in general are assumed to be anisotropic, with pa- 
rameters (3\ and The spins are assumed to be clas- 
sical variables and are constrained to lay on the plane. 
Measurements taken in a considerable variety of layered 
compounds reveal that the inplane spin exchange inter- 
actions are isotropic to a high accuracy, meaning that the 
magnetic crystalline anisotropy dominates- 1 ' 2 . Therefore, 
we consider (3\ ~ f3 2 and adopt a% and a 2 as the source 
of anisotropy in our model. 

To take the continuous limit of Eq. we partition 
the antiferromagnetic lattice in two ferromagnetic sub- 
lattices, and define a mean field for each sublattice as 
m||(r) = (riiS^g^ , calculating the spin average in a 
ball centered at the point r, in real space. In terms of 
the staggered fields (m-f, m^), we define a continuous spin 
OP as 

m(r) = (m T - m x ) (r) , 

in the limit when the lattice constant goes to zero. Its 
norm m is a positive number between and 1. In the 
absence of dynamics, ixi| and nij are locally antiparallel 
and thus m is maximum in the saturation regime, for 
T <C Tn ■ We also define a renormalized spin occupation 
density 77, varying in the interval [0, 1], which plays the 
role of the secondary OP (it does not vanish at Tjv). To 
remove the criticality at T/v, we impose 



where t = 1 — T/Tn is the reduced temperature, and 
f3 > corresponds to the OP critical exponent. If we con- 
sider now the OP components along the directions par- 
allel and transverse to the stripes, meaning m\\ = mco$9 
and m± = msm.9, and use relation JSJ, we find that the 
macroscopic regime is well defined in terms of the two 
OP degrees of freedom, 77 and 9. The symmetry group 
of the OP is the subgroup C 2v of 0(2), which includes 
inversion of the spins and two reflection planes along the 
two principal axes of the magnetization. The continuous 
limit of the microscopic Hamiltonian J2J leads to the GL 
free energy density (see the Appendix) 



(t 2 ?b -a )(v) 2 
[a 2 + c 2 t 213 - t 2p b 2 (ay cos 



t 2 Pb {8) 2 r] 2 - 2a 2V 



a± sur 



(4) 



The dot means a spatial derivative with respect to z. We 
assume a 2 and c 2 to be regular functions of temperature, 
say a 2 (T) and c 2 (T), while b 2 = bt is a standard Landau 
second order term parameter, with b a positive constant. 

To satisfy the continuity condition of the magnetic field 
in the stripe interface with the antiferromagnetic environ- 
ment, we must have m-ej^ = inside the bulk of domains. 
It leads to fixing the direction of easiest magnetization 
along §|| , what is just equivalent to impose a± < an to 
the anisotropy parameters. The C 2v point group of the 
spin ordered phase admits two invariant terms, m 2 and 
m 2 cos(2#), which generate the polynomial expansion of 
the free energy potential^. The second order anisotropy 
term b 2 m 2 (a± sin 2 9 + at\\ cos 2 9) can be decomposed in 
such basis as ifo 2 {(an + ckjJto 2 + {a\\ — a±)m 2 cos(29)} , 
making evident the competition between the full symme- 
try term of the paramagnetic phase (0(2) point group) 
and the C 2v symmetry one. From this, the proper ex- 
pression for the inplane crystalline anisotropy is 



< e 



< 1 



(5) 



The appropriate set of boundary conditions (BC) that 
applies to the free energy (0J is 



0. 



z — > —00 

7T Z — > OO 

Z — > ±00 



and 



r\ = 1, z 
fj = 0, z 



±00 
±00 



(6) 
(7) 

(8) 
(9) 



The SDW modulates transversely to the stripes, in or- 
der to satisfy p = 1 — r\. If we denote the energy of a 
single domain by E = f dz , variation under the 
constrain of the global charge conservation QJ, yields 



Tj(r) = m(r) t 



(3) 



S(E + AU) = 6 d 



r],V,r],z) 



A 77 = 



4 



where A is a Lagrange multiplier. The variational equa- 
tion above results in 



r)L ~ lj where r\h satisfies the equation 1)14(1. we en- 
counter that 



tanh (A z) 



(10) 



1>(z)~g(z) [sech(A-^)] 



(16) 



and 



Aij - (5- C/sech^A- 1 z)) T) - 3C(1 - rf) + 

+2L>(l-?7 3 ) + 5 = , (11) 

where A -1 = [(an — a±)b 2 /bo]^ is the domain wall width 
in the spiral order. Similarly to Bloch walls in ferro- 
magnets, a finite A results from the competition between 
the exchange and the anisotropy terms in the magnetic 
energy^. In expression (|11[> . we have the quantities 
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u 
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D 



b t 2P - a , 

a 2 — i 2/3 6 2 «|| +c 2 t 213 

-2{a\\-a 1 _)b 2 t 2li , 



where g(z) is a smooth and convergent hypergeometric 
function 

g(z) = f(a_,A+; 1 + xW; sech 2 (A _1 z)) , 
with A± = r± + \\fV and 

V = ^x[2D(3- x)-3C] 

as an explicitly doping dependent quantity. 

Going one more step further, we apply the constraint 
JTJ to the perturbed solution t}l + ci/j, which yields 



J_>[1-^l( U )]-^ 
J™ duip(u) 



(17) 



c 2 i 2/3 + t* p b. 



which are defined in terms of the original GL parameters 
given in J3J. As usual, ciq and bo have no dependence 
on temperature and are assumed to be independent of 
(o,2, b 2 , c 2 t 2 P ,b4t 4 @) in the parameter space. The limit 
ao,&o — > 0, with A -1 — > oo, reproduces the Ising case 
where (|ll|l becomes a z-independent equation 



where u = z/X is the coordinate scaled by the soliton 
width. 



A. Estimation of parameters 

A numerical calculation shows that J_ du (1 — — 
0.81 UX 2 /A, with u — z/X. In lowest order we write the 



5(1 - 77) - 3C(1 - T] 2 ) + 2D (1 - t? 3 ) = , (12) constraint (TJ as /du (1 - 77^) = fx, and find that 



which yields a constant r\ = 1 — x. Substituting 77 in 
terms of the doping x in i|12[) one finds 



d 4(7 „ 
A 5 A 



b t 



5 do — bot 



(18) 



5 = 3C (2 - z) - 2/J [3(1-. 



(13) 



for x 7^ 0. In order <|13|l to be valid in a continuous inter- 
val of temperature, we must have a 2 (T) cx t 2 , c 2 (T) cx t 
and the OP critical exponent f3 = i. 

Now we solve equation <|11[) in a perturbative way, in 
relation to the parameters (c 2 ,bi) which are associated 
to the cubic and quartic terms in the free energy density 
(gj. We write 

v( z ) = VL{z) + cip(z), 

with tjl satisfying eq. (|11|) in lowest order and if> cor- 
responding to higher order corrections due to c 2 and 64, 
which are small parameters in the GL sense. Imposing 
the limit c 2 , &4 — > 0, condition 1)1 3[) becomes 5 = for 
x 7^ 0, and we obtain from ifTTl) 

AfjL + J7sech 2 (A _1 z)t]l = , (14) 

whose solution is the hypergeometric function 

77L(z)=^(L_,r + ; l;sech 2 (A- 1 z)) , (15) 



Up to this order, the amplitude of the CDW peaks de- 
pends on the adimensional parameter UX 2 /A only. This 
latter quantity is numerically limited to the interval [0, 2] 
in order to have < p(z) < 1. As a consequence, relation 
(|18|l leads to the approximate inequality 

„ d / 8 

The above condition also allows us to estimate the limi- 
tation of the temperature to saturate the peak. Within 
our simplified model we get ti ow ~ (ao/26o) < 1, when 
the peaks have maximum amplitude (for UX 2 /A ~ 2). 
Two examples are shown in Fig. 2. 

We observe that the limit d/X — ► 1+ for the ratio be- 
tween the domain width and the domain wall width, is 
indicative of a CDW commcnsurability transition from 
the present low incommensurable (LI) phase to some in- 
commensurable one. Because the validity of this LI phase 
model requires that d > A, we get the condition 



x < 



b t 



5 cio — bot 



(19) 



with r± = |(1 ± y/l + 4UX 2 /A). Solving the linearized 
V> equation resulting from (|ll|l in the asymptotic limit 



which delimits region LI from above in the phase diagram 
displayed in Fig. ^ 
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FIG. 2: CDW solution for a /&o = 1.5. On the left: U\ 2 /A = 
1 at t = 0.58, with arbitrary parameters xd/X = 1 and FA 2 = 
0.2. We define an arbitrary scale Ao for the peak width at this 
temperature. On the right: rescaling the reduced temperature 
to t = 0.31, we have U\ 2 /A = 0.25, xd/X = 0.25 and VX 2 = 
0.05, with A = 1.69 A . 



III. STRIPE PHASE TRANSITION 

In the spirit of GL commensurability transitions, we 
now pursue the phenomenological description of the 
stripe critical regime, where magnetic order is absent. 
For high incommensurability, the appropriated OP is 
Spe 1 ^, which describes the charge modulation around the 
uniform high- symmetry phase charge distribution pq. In 
the vicinity of the critical point Tg, the modulation is 
soft and therefore we assume that (Sp) = 0, and we de- 
scribe the CDW as a periodic modulation with constant 
amplitude p(z) — po = Spe^^. 

The generalized GL free energy density is written as a 
contribution of an explicitly dependent term on the OP 
derivatives K plus a potential U(tp). The most general 
"kinetic" term has the form 



K 



(20) 



The invariant part of the antisymmetric term in i|20|) 
is called the Lifshitz invariant, and is the essential in- 
gredient to stabilize a modulated incommensurate OP 
phase 13 . The Lifshitz invariant comes from the most gen- 
eral expansion which includes first order derivatives, af- 
ter one removes contributions which are not extensive in 
character (terms which contribute to the total free energy 
F as surface integrals). 

The point group of the OP is chosen with a low sym- 
metry, just noting that a sinusoidal charge distribution 
locally displays mirror symmetry in the phase <f) at max- 
ima and minima of the distribution. Thus, with an ap- 
propriate phase choice, we include in the symmetry group 
a reflection a to describe the operation <j> — * — (f) and the 
point group is C± v = {e,c}, where e is the identity^*. 
The real components of the OP, (Sp cos 4>, 5psm<fi), trans- 
form differently: Sp cos <f> is invariant and Spsinfj) trans- 



form after the antisymmetric representation. Therefore, 
the polynomial basis to generate invariants of the group 
is {(dp) 2 , dp cos (j)}. We observed that the latter basis 
displays particle-hole symmetry, since 

—Sp cos 4> = Sp cos (4> — tt) = Sp cos <f>' , 

i.e. the change in sign of the charge may be included in a 
redefinition of the phase variable (interchanges maxima 
and minima of the charge distribution, leaving the free 
energy invariant). 

If we write the 'kinetic' energy (|20|l using the real com- 
ponents of the OP, and using the fact that (Sp) — 0, we 
get K in its most symmetric form, with functions in front 
of the derivatives which are invariants of the symmetry 
group. Note that (Sp) 2 cos 2 <j) is not included, if we addi- 
tionally require the 'kinetic' term to be invariant under 
an arbitrary phase shift: 



-a O i(5p) 2 q> + a o (Sp) 2 (0) 2 



(21) 



The above 'kinetic' Hamiltonian l|21|) represents the elas- 
tic energy due to the distortion of the OP in relation to 
the homogeneous case. The phenomenological parameter 
doi is related to the Lifshitz invariant, which is linear in 
(f>. It is then clear the role of this term: while 'kinetic' 
and 'potential' energies compete, the Lifshitz invariant 
will try to stabilize a phase with a non vanishing <fi. 

The 'potential energy' is expanded as a power series of 
the polynomial basis of invariants {(Sp) 2 , Sp cos (/)}. This 
way, the most general form of the potential up to third 
order is 

U = a±8p cos0 + a 2 (5p) 2 +a 2 i(5p) 2 cos 2 + 

a 3 (Spf cos0+ a 31 (5p) 3 cos 3 <j) + a^Spf (22) 

In first place, note that particle-hole symmetry (which 
has been used in the choice of the base) , excludes a term 
of the type (Sp) 3 cos 2 <f>. This is a key point in relation 
to the order of the transition at Tg. More on this later. 
As usual in GL theories, the second order coefficients a 2 
and a 2 \ are assumed to be strongly temperature depen- 
dent. At last, concerning the fourth order terms in 122|l . 
a comment is pertinent. Our aim in this section is to 
describe the transition between the harmonic incommen- 
surate phase (HI in Fig. 1) and the high symmetry phase 
(HS in Fig. 1), where the OP amplitude Sp is small (note 
that the Lifshitz invariant term is proportional to (Sp) 2 ). 
A single fourth order term is then included to provide 
for a global stability of our solution, assuming 0,4 > 0. 
Note that a different physics may be incorporated with 
other fourth order terms. Take for instance the case of 
(Sp) 4 cos 4 </>, with a negative coefficient. Clearly, this is 
a 'lock-in' term that should favor a homogeneous com- 
mensurate solution for large Sp. For simplicity, this case 
will not be considered here. 

The free energy (K + U) does not explicitly depend 
on the variable z. In this case, the variational Euler- 
Lagrange minimization leads to a first integral of the form 
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A = (K + U) - <j> 



d(K + U) 



(23) 



where A is a constant. The explicit calculation of i|23|) 
results in the equation 



= n 



a 21 cQg 2 ( 



which couples <j) with Sp, with the definitions 
v _ + a 3 (Jp) 2 



and 



17 = 



a 5p 



ai5p + a 2 {5p) 2 + a 3 (Sp) 3 + a 4 (Sp) 4 
a (Sp) 2 



A 



(24) 



(25) 



(26) 



which are constants that parametrize our solutions. 

If we take the limit 021, 031 — * in l|24l) . the remaining 
expression is a time-independent sine-Gordon equation, 
which describes highly incommensurable solutions in the 
limit v/VL <C 1 (note that v/Vl ~ Sp). For finite but small 
corrections in (021,031), we get 



0/2 



dtp 



a-21 



1 



4a r2 



3 u 
4ft 
3 



(f>+ - sin(2</>) 



a 3 i5p 



(27) 



The left hand side of (|27|l is an elliptical integral of the 
first kindi^. For 021, 031 — > 0, </>(z) = 2am(u) ~ 
fi 1 / 2 ^ = fc 2, for |j -C 1, where am is the amplitude 
of the Jacobian Elliptic function. Carrying on this pro- 
cedure perturbatively up to terms of the order of (Sp) 3 , 
we get 

cf>(z) ~ kz + ( 1 + ~\ sm(2k z) + 



8a a n V 4 fi 
1 



kz + e(z) 



x ( 1 - -snT (k z) 



(28) 



with 



fc ~ kt] 



a-21 



Aa Q \ + 4 ft 



From this we can compute the CDW energy E — 
I J Q Z dz (K + U) and obtain k and Sp which stabi- 
lize the stripe phase at the global minimum condition 
dE/d(Sp) = dE/dk = 0. Since ft ~ 1/ip 2 , the terms 
a 2 i/fl and a^iSp/fl are of orders 2 and 3 in Sp, respec- 
tively. This implies that the only term that carries cor- 
rections from e of the order of (Sp) 3 in the free energy, 
is a\Sp cos (fcz + e) , which expanded within this order 
becomes 

a\5p [cos(fcz) — sin(/cz) e] = a\5p cos(k^z) + 



01 



sin(koz)sin(2k z) + o(Sp' i ) (29) 



For l z ^> 2ir/ko, the integrals f "dz {sin(fcoz) sin(2koz)} 

and f "dz cos n (kz), with n odd, are non-extensive quan- 
tities in the energy and may be ignored. The energy will 
be given by the simple expression 



E = IL 



(M 2 (- 



a ik + a k + a 2 



a-21 
2 



a 4 (Sp) 4 



where 04 is a renormalized fourth order coefficient, which 
includes small contributions from 031 and 021. Therefore, 
we assume that the condition 04 > is fulfilled. Note 
that third order terms in Sp do not contribute to E in the 
thermodynamic limit. The expression l|30() is a minimum 
for 



goi 
2a 



and 



Sp = 



qgi - 2ao(a 2 i + 2a 2 ) 



(31) 



(32) 



In the limit of highly incommensurable solutions, the 
potential U is irrelevant for the stabilization of the 
ground state. The stability is reached through the com- 
petition of the ao squared gradient term (which usually 
fixes the scale of energy in such classes of phenomenolog- 
ical models), with the Lifshitz invariant coefficient doi, 
what we interpret as a minimization of the CDW elastic 
energy. 

When T — > T$ from below, Sp — * 0, and the Lifshitz 
invariant coefficient aoi competes with the second order 
coefficients (a 2 , a 2 \) of the potential U , meaning they 'co- 
operate' with ao. But far below the stripe critical point 
T<j, the importance of the Lifshitz invariant term is small 
in comparison to the 'potential' energy contribution of 
the several cosine terms, which are known to stabilize LI 
solutions, as we have shown in Section II. Therefore, the 
second order coefficients (02,021) should change sign in 
between of these two limits. Thus, the precedent analysis 
requires the coefficients a 2 and 021 to be defined in terms 
of a new critical parameter r = T/Tq — 1,Tq<Ts, which 
drives a CDW commensurability transition. 

Writing a 2 — ar and 021 = o! r (a, a! > 0), the stripe 
phase transition at T5 will be of second order if Tq and 
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T$ are related by the constraint 



T s = T 1 



2a (a' + 2a) 
Replacing (|33|) in i|32|) we find 



(33) 



'2a (a' + 2a)+c4 (T a -T 



(34) 



which goes to zero at T$ under the usual mean field 
square-root law. 

The periodic e(z) function in l|28|l introduces non- 
periodical corrections in the CDW phase (f> — kz + e(z). 
The only way to preserve the CDW periodic nature is 
to set e(z) = 0, i.e et2i, Q.31 = 0. In this picture, the 

CDW solution is just p(z) = x + 5pcos(- — z), using ap- 



propriate units, with x being the doping. Applying the 
global constraint I L" p(z)dz = Q in the limit l z ^ 3?, 

we encounter that po = x. Since Sp < po, the (Sp) ~ 
region is delimited by the inequality x > Sp in the phase 
diagram depicted qualitatively in Fig. 1. 



since there are no magnetic correlations to compete with 
the repulsion of holes. For this phase, we construct the 
free energy from general theoretical arguments coming 
entirely from symmetry. We conjecture that this free 
energy includes, in a stationary fashion, dynamical phe- 
nomena that occur at high temperature, such as elastic 
and electronic coupled effects that induce a polarization 
of the lattice. In an actual phase diagram of cuprates, as 
the one shown in Fig. 4 of Refll7l this phenomenon is 
located in the high temperature sector, considerably far 
from the tetragonal-orthorhombic structural phase tran- 
sition. This latter is still in the interior of the 2D antifer- 
romagnetic correlated region, and has no relation with 
the high temperature physics that we are dealing with 
here. 

The C\ v point group of the OP for the HI phase implies 
the occurrence of odd terms in the free energy density. 
However, particle-hole symmetry makes those terms non- 
extensive in character, and they can be neglected in the 
thermodynamic limit. This means, in principle, that one 
may not rule out the possibility of a second order phase 
transition for stripes, in contrast to other phenomenolog- 
ical approaches^. 



IV. DISCUSSION 

We present here a brief summary of the paper and we 
discuss the results of the last two sections. Phase I de- 
scribes the LI spin driven regime of stripes, which takes 
place in our phase diagram for low doping concentra- 
tions. This phase is also experimentally associated with 
incompressible stripes (for x < |, in the cuprate case). 
This is an evidence that they interact weakly, despite 
the existence of 'long-range' spin-spin and charge-charge 
interactions (in both cases, the correlation length is of 
the order of few hundred angstroms). Separately, each 
of them would drive the system to macroscopically ho- 
mogeneous phases. However, their competition, even in 
the absence of dynamical correlations, reduces consider- 
ably the effective range of the interactions and stabilizes 
non-homogeneous phases in a smaller scale. 

We showed in Section I that a simple static Hamilto- 
nian, like the one given by ||2J), captures in a qualita- 
tive way the physics of competing effects which leads to 
low temperature CDW formation in the presence of two- 
dimensional magnetism. The second conclusion is that 
this soliton-like CDW scenario changes drastically with 
temperature. Moving in region LI at constant doping 
and increasing temperatures (see Fig. 1), our solution 
shows that the ratio between the inter-stripe distance 
and the stripe width decreases (as well as the amplitude 
of the peaks) and goes to the limit d/X — > 1, which cor- 
responds to the commensurability transition that sepa- 
rates the LI phase from the high incommensurable (but 
still spin driven) phase I. 

In contrast, in the absence of 2D antiferromagnetism, 
we may ask which is the mechanism behind the CDWs, 
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Appendix 

In this appendix we will demonstrate that the GL free 
energy is equivalent to Hamiltonian (J2J) in mean field 
approximation. Starting from the definition of our stag- 
gered mean field m, we calculate macroscopic average of 
the spin term in @, which yields 

2 

H s = - d 2 r ajm^r) 
2 

- f { d 2 rd 2 e J'J2 / WW* + £ ) > (35) 

in the limit where the lattice parameter goes to zero. The 
charge term calculation is similar, and we may write 

H c = //d 2 rd a e ± [1 - V (r)} [1 - r,(r + e)} 

Observing that the OP is by definition a continuous func- 
tion that varies slowly in space, we assume that r](r + e) 
and m(r + e) can be expanded up to terms of lowest or- 
der. Performing the integration in e, the resulting free 
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energy is a functional of the form 

2 

/(r) = -a (r)) 2 + 02(1 - r?) 2 + ^ [M™j) 2 - b 2 a 3 m 2 ] , 

(36) 

with oij taken as an effective anisotropy parameter. 

Note that the charge (1 — 77) acts as a secondary order 
parameter, with no critical behavior. In order to achieve 
a minimum for the thermodynamic potential, we must 
include higher order powers of the primary OP m in the 
above expression. Since the odd ones break the free en- 



ergy invariance under reflections, the lowest correction is 
of fourth order. So, we heuristically introduce the sim- 
plest term 64 m 4 , which represents an isotropic magnetic 
quadrupole interaction. In the same way, we may con- 
sider a coupling term for the spin and charge degrees of 
freedom. We rule out the OP coupling with odd powers 
of (1 — 77), due to particle- hole symmetry. If we assume 
that the coupling does not depend on the crystal direc- 
tions, the lowest order term is ci (1 — rj) 2 m 2 . This closes 
expression (0J in our model. 
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